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SPECIFICATION 



Radio-Wave Arrival-Direction Estimating Apparatus and Directional 

Variable Transceiver 

FIELD OF THE INVENTION 

The present invention relates to a radio-wave arrival-direction estimating 
apparatus employing an array antenna, and a directivity variable transceiver for 
varying antenna directivity based on an estimation result from the estimating 
apparatus. 

BACKGROUND OF THE INVENTION 

An arrival direction of radio wave is conventionally estimated accurately in 
a method such as Multiple Signal Classification (MUSIC) method, using an array 
antenna comprising a plurality of antenna elements. The MUSIC method is 
disposed in R. O. Schmidt, "Multiple Emitter Location and Signal Parameter 
Estimation", Institute of Electrical and Electronics Engineers (IEEE) Trans., AP- 
34, pp. 276-280 (1986). This method includes an algorism for accurately 
estimating a direction of a plurality of incident waves with the same frequency 
band. 

In this method, M (>1) antenna elements receive signals, and a receiving 
unit connected to each antenna element converts the frequency of each of the 
received signals, detects a phase of it, and demodulates the received signal to a 
signal comprising orthogonal I and Q signals. An analog/digital converter (A/D 
converter) converts the demodulated signal to sampling data and outputs the 
data to a direction estimating processor. The direction estimating processor 



estimates a direction of the incident waves using the sampling data by the 
MUSIC method. In other words, using sampling data xi (k), X2 (k), ... , xm (k) at 
sampling time kAT obtained by respective antenna elements, a correlation matrix 
calculation unit creates receiving vector x (k) written as 
5 x(k)-[x x (k)x 2 (k)...x M (k)J (Equation 1), 

where T shows transposition of a vector, AT is a sampling interval, and k is a 
natural number. The correlation matrix calculation unit, using receiving vectors 
x (k) for k = 1 to N, further finds MxM correlation matrix R written as 

R = ^x(k)x(k) H IN (Equation 2), 

10 where H shows complex conjugate transposition of a vector. 

The calculation unit calculates eigenvalues K\ - Xm of correlation matrix R 
in the descending order, and eigenvactors ei - eM corresponding to eigenvalues X\ - 

Next, the calculation unit calculates an evaluation value of an arrival-angle 
15 evaluation function, assuming number of the incident waves is S, and using noise 
spatial eigenmatrix En = [es+i, ... , eM] and a feature that signal eigenvector space 
Es = [ei, ... , es] and En are orthogonal to each other. This En is formed with (M 
- S) eigenvactors, namely column vectors, belonging to a noise eigenvactor space 
having the relation written as 
20 A 2 ^"-^A 5 >A s + i = A s + A M (Equation 3), 

and Es is formed with eigenvactors ei - es. In other words, arrival-angle 
evaluation function F(8) for evaluating orthogonality between En and Es is 
defined by 

" —h (Equation 4), 

25 where a(6) is a complex response (hereinafter called a steering vector) of the array 
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antenna as a function of azimuth 0. Azimuth 0 varies in a predetermined angle 
range. When azimuth 0 equals to the arrival angle, ideally, arrival-angle 
evaluation function F(0) is infinite. A resultant peak direction of F(0) from 
calculation for the varied 0 is set to be the arrival-angle evaluation value of the 
5 incident waves. 

Number S of incident waves is generally unknown, so that the number is 
determined based on an eigenvalue distribution and number-of-signal 
determination criteria. The criteria is described in M. Wax and T. Kailath, 
"Detection of Signals by Information Theoretic Criteria", IEEE Trans. On 
10 Acoustics, Speech and Signal Processing, Vol. ASSP 33 (2), pp. 387 - 392, 
February (1985). 

The radio-wave arrival-direction estimating apparatus employing the 
MUSIC method discussed above estimates an arrival direction accurately by 
signal processing, using an algorithm of deriving engenvalue of a correlation 

15 matrix of array received signals; In such an estimating apparatus, correlation 
between waves generated by reflection on the ground or a building increases 
when a relative delay time between these waves is shorter than a symbol length. 
In this case, correlation matrix R degrades, and therefore the incident waves 
cannot be precisely separated. 

20 For preventing the degradation, a spatial smoothing technique is proposed. 

This spatial smoothing technique is described in Pillai et al, "Forward/Backward 
Spatial Smoothing Techniques for Coherent Signal Identification", IEEE Trans. 
On Acoustics, Speech and Signal Processing, Vol. 37, No. 1, 1989. The example 
has estimated the arrival direction using spatial samples from the array antenna; 

25 however the MUSIC method can be similarly applied to a signal sampled every 
frequency and the delay time of the received waves can be estimated at high 
resolution. 
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The estimation accuracy of the arrival direction in the MUSIC method 
depends on variation step &6 of 6 in the arrival-angle evaluation function (Eq.4). 
When A0 increases, a calculation amount in the entire variation range of 9 
decreases, but the peak direction of the arrival-angle evaluation function cannot 
5 accurately detect the peak direction and the accuracy decreases. When A9 
decreases, the peak direction of the arrival-angle evaluation function can be 
accurately detected, but a calculation amount in the entire variation range of 8 
increases. 

10 SUMMARY OF THE INVENTION 

It is an object of the present invention to provide a radio-wave arrival- 
direction estimating apparatus allowing reduction of a total calculation amount 
for an arrival-angle evaluation function without causing accuracy degradation of 
the arrival direction. It is another object of the present invention to provide a 

15 directivity variable transceiver for improving transmitting and receiving quality 
by controlling antenna directivity. 

In the present invention, product of a noise spatial eigenmatrix and a 
conjugated and transposed matrix of it is a product of an upper or lower 
triangular matrix. Therefore, the calculation amount for the arrival-angle 

20 evaluation function can be reduced in the entire angle range for the estimation of 
the arrival direction. The arrival angle evaluation using the arrival-angle 
evaluation function that has a heavy calculation load can be significantly reduced 
during angle sweeping in the MUSIC method. Processing of the arrival- 
direction estimating apparatus can be speeded or an apparatus structure can be 

25 simplified. 
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Fig. 1 is a block diagram of a radio-wave arrival-direction estimating 
apparatus in accordance with exemplary embodiment 1 of the present invention. 

Fig. 2 is a graph illustrating reduction of a calculation amount in 
accordance with exemplary embodiment 1. 
5 Fig. 3 is a graph illustrating a calculation amount required for cholesky 

factorization in accordance with exemplary embodiment 1. 

Fig. 4 is a block diagram of a direction estimating processor in accordance 
with exemplary embodiment 2 of the present invention. 

Fig. 5 is a block diagram of a direction estimating processor in accordance 
10 with exemplary embodiment 3 of the present invention. 

Fig. 6 is a block diagram of a direction estimating processor in accordance 
with exemplary embodiment 4 of the present invention. 

Fig. 7 is a block diagram of a direction estimating processor in accordance 
with exemplary embodiment 5 of the present invention. 
15 Fig. 8 is a block diagram of a direction estimating processor in accordance 

with exemplary embodiment 6 of the present invention. 

Fig. 9 is a block diagram of a direction estimating processor in accordance 
with exemplary embodiment 7 of the present invention. 

Fig. 10 is a block diagram of a direction estimating processor in accordance 
20 with exemplary embodiment 8 of the present invention. 

Fig. 11 is a block diagram of a direction estimating processor in accordance 
with exemplary embodiment 9 of the present invention. 

Fig. 12 is a block diagram of a direction estimating processor in accordance 
with exemplary embodiment 10 of the present invention. 
25 Fig. 13 is a block diagram of a direction estimating processor in accordance 

with exemplary embodiment 11 of the present invention. 

Fig. 14 is a graph illustrating an operation of a direction estimating process 



6 

in accordance with exemplaxy embodiment 11. 

Fig. 15 is a block diagram of a directivity variable receiver in accordance 
with exemplary embodiment 12 of the present invention. 

Fig. 16 is a block diagram of a directivity variable transceiver in 
5 accordance with exemplary embodiment 12. 

Fig. 17 is a block diagram of a directivity variable receiver in accordance 
with exemplary embodiment 13 of the present invention. 

Fig. 18 is a block diagram of a directivity variable . .transmitter in 
accordance with exemplary embodiment 14 of the present invention. 

10 

DESCRIPTION OF THE EXEMPLAHY EMBODIMENTS 
Exemplary embodiments of the present invention are demonstrated 
hereinafter with reference to the accompanying drawings. 

15 1. FIRST EXEMPLARY EMBODIMENT 

Fig. 1 is a block diagram of a radio-wave arrival- direction estimating 
apparatus in accordance with exemplary embodiment 1 of the present invention. 
Array antenna 1 comprises M (>1) antenna elements 1-1 to 1-M. Antenna 
elements 1-1 to 1-M receive high frequency signals 2-1 to 2-M. Receiving units 

20 3-1 to 3-M connected to antenna elements 1-1 to 1-M convert frequency of the 
signals and demodulate the converted signals to signals 4-1 to 4-M comprising 
orthogonal I signal and Q signal. AID converters 5-1 to 5-M sample respective I 
signals and Q signals of the demodulated signals 4-1 to 4-M, and convert the 
demodulated signals to complex digital signals 6-1 to 6-M. Each of the complex 

25 digital signals has the I signal in its real part and the Q signal in its imaginary 
part. 

Correlation matrix calculation unit 7 creates receiving vector x (k) given by 
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Eq.l, using complex digital signals xi (k), X2 (k), ... , xm (k) at sampling time kAT 
derived from complex digital signals 6-1 to 6-M. Here k is a natural number and 
AT is a sampling interval. Correlation matrix calculation unit 7 further derives 
M x M correlation matrix R written by Eq.2 using receiving vectors x (k) 
5 accumulated for N sampling periods. 

Noise spatial eigenmatrix calculation unit 8 applies eigenvalue 
factorization to derived correlation matrix R to derive eigenvalues k\ - Xm in the 
descending order and eigenvactors ei - eu corresponding to them. When a 
number of the incident waves is S, calculation unit 8 outputs noise spatial 
10 eigenmatrix En = [es+i, , eu] comprising (M - S) eigenvactors, namely column 
vectors, belonging to a noise partial space having the relation given by Eq.3. 

Triangular matrix calculation unit 9 derives product U of a noise spatial 
eigenmatrix and a conjugated and transposed matrix of it as in, 

U = En En h (Equation 5). 

15 Because matrix U is an M x M positive definite matrix, triangular matrix 

calculation unit 9, using cholesky factorization, derives lower triangular matrix L 
written as 

U=LL H (Equation 6). 

Arrival-angle evaluation unit 10 evaluates an arrival angle every 
20 predetermined angle step A0 using the evaluation function 

Fl{6) = a H (d)LL"a(8) _ 

^ W (Equation?), 

= \\L H a{B) f 

where, | I x | | is the norm of vector x, and a(0) is a normalized steering vector of 
the array antenna. This evaluation function (Eq.7) is derived from the arrival- 
angle evaluation function written as Eq.4 using lower triangular matrix L. 
25 Since elements outside the lower triangular part in lower triangular matrix L are 
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null, a ratio of a sum-of-product calculation amount for the arrival-angle 
evaluation function given by Eq.7 to that in Eq.4 is (M+3) / [2 (M-S+l)]. The 
calculation amount for the function given by Eq.7 can be therefore reduced if 
number S of incident waves satisfies S < (M-l) / 2. 
5 Referring now to Fig. 2 there is shown a ratio of the calculation amount for 

the arrival-angle evaluation function given by Eq.7 to that by Eq.4 when number 
S is .1. Fig. 2 shows that a reduced calculation amount in Eq.7 in this method 
increases as the number of antenna elements increases. For example, when the 
number of antenna elements is 6, the calculation amount is about 75% of that in 
10 the prior art. 

Referring now to Fig. 3 there is shown a ratio of the calculation amount for 
the cholesky factorization to the calculation amount for arrival-angle evaluation 
function F(9i) for one arrival angle 9i in Eq.4 . Fig. 3 shows that the calculation 
amount for the cholesky factorization does not reach the calculation amount for 

15 arrival-angle evaluation function for five arrival angles even if the number of 
antenna elements is 20. The arrival angle evaluation is usually performed for 
more than 5 arrival angles, so that the increase of the calculation amount for the 
cholesky factorization can be considered extremely smaller than that for the 
arrival angle evaluation in the entire range of the angle sweep in Eq.7. 

20 Arrival-angle determination unit 11 detects a peak direction based on an 

arrival-angle evaluation result every AG in a variable range of 8, and uses the 
detected direction as an arrival-angle estimation value of the incident waves. 

In the present embodiment, triangular matrix calculation unit 9 applies 
the cholesky factorization to product U of the noise spatial eigenmatrix and the 

25 conjugated and transposed matrix of it to derive lower triangular matrix L. 
Arrival- angle evaluation function F2(8) is derived by equivalent conversion of the 
arrival-angle evaluation function written as Eq.4 using matrix L. Additionally, 
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using function F2(8), the calculation amount for the arrival angle evaluation can 
be reduced under the condition S < (M-l) / 2. 

The direction estimation using the MUSIC method has been discussed 
above. However, the method of the present embodiment can be as-is applied to a 
5 received signal sampled eveiy frequency when the delay time of the received 
wave is estimated at high resolution. That is because an evaluation function 
similar to that in Eq.4 is used. 



Eq.6; however, an upper triangular matrix also obviously produces a similar 
10 advantage. 

Triangular matrix calculation unit 9, using a modified cholesky 
factorization, may find lower triangular matrix L and diagonal matrix D 
expressed as 



The modified cholesky factorization does not require the square root calculation, 
so that calculation time can be reduced. 

In this case, arrival-angle evaluation function F<2(6) is expressed as 



The lower triangular matrix is derived by the cholesky factorization in 



U=LD L H 



(Equation 8). 



15 



F a (0)- 




(Equation 9), 



20 where bk and dk are vector elements given by 



(Equation 10) 




and 
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Id 




0 \ 



D = 



0 



(Equation 11), 




0 d 



respectively. The lower triangular matrix has been used in the present 
embodiment; however, an upper triangular matrix also obviously produces a 
similar advantage. 

Additionally, correlation matrix calculation unit 7 may apply a spatial 



wave. The method in the present embodiment can be similarly applied in this 
case if the spatially smoothed correlation matrix instead of correlation matrix R is 
fed into direction estimating processor 12. 

An example of an array antenna having a constant-interval linear' array 
shape is described in M. Haardt and J. A. Nossek, "Unitary ESPRIT: How' to 
Obtain Increased Estimation Accuracy with a Reduced Commutational Burden/' 
IEEE Trans. Signal Processing, vol. 43, No. 5, pp. 1232 - 1242 (1995). In this 
example, thanks to the conjugation center symmetiy of a phase of a steering 
vector, the steering vector can be converted into a real vector using unitary 
matrix Qm written as 



where a(B) is a steering vector when the phase center matches with the array 
center. 

A method similar to the method in the present embodiment can be applied 
to this case, if direction estimating processor 12 uses the real part of Qm h RQm 
instead of correlation matrix R, and b(0) instead of steering vector a(0). 

When an array antenna having a linear array shape is employed, 
estimation accuracy in the end fire direction decreases. Therefore, arrival-angle 
evaluation unit 10 sets the angle interval in the end fire direction of the array 



smoothing technique to the correlation matrix in order to suppress correlation 



b{6)-Q H M a{e) 



(Equation 12), 
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antenna to be larger than that in the bore-sight direction, and calculates an 
evaluation value of the arrival-angle evaluation function. The calculation 
amount can be thus reduced. Here the bore-sight direction means the direction 
of the normal to the linear array arrangement direction. 
5 A configuration may be employed that comprises, instead of receiving units 

3 and A/D converters 5 in the present embodiment, the following components: 

intermediate frequency (IF) receiving units for converting frequencies 
and detecting phases of RF signals supplied from respective antenna 

elements 1-1 to 1-M and outputting IF signals; 
10 IF A/D converters for converting the IF signals to digital signals; and 

digital orthogonal wave detectors for orthogonally demodtriating the 
digital signals and supplying complex digital signals to the correlation matrix 
calculation unit. 



15 2. SECOND EXEMPLARY EMBODIMENT 

Fig. 4 is a block diagram illustrating the other configuration of direction 

estimating processor 12 of the radio-wave arrival- direction estimating apparatus 

of the present' invention. Components other than the direction estimating 

processor in the radio-wave arrival-direction estimating apparatus are similar to 
20 those in embodiment 1 described in Fig. 1, so that diagrams and descriptions of 

these components are omitted. Components different from embodiment 1 will be 

mainly described hereinafter. 

Operations until correlation matrix R is fed into direction estimating 

processor 12a are similar to those in embodiment 1. Inverse matrix calculation 
25 unit 20 calculates inverse matrix R 1 of correlation matrix R. Because inverse 

matrix R l is a positive definite matrix, triangular matrix calculation unit 21 

derives lower triangular matrix L written as 
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R' l =LL H (Equation 13), 

using the cholesky factorization. 

Using lower triangular matrix L, arrival-angle evaluation unit 22 converts 
a conventional arrival-angle evaluation function obtained by Capon method as 
5 given by 

Fc( 6 ) . — — (Equation 14), 

a" \d)R- l a{6) h 

to an evaluation function 

FCl(e) -jiI%W (Equation 15), 

where, I I x | | is the norm of vector x, and a(6) is a normalized steering vector of 
10 the array antenna. Arrival-angle evaluation unit 22 then evaluates an arrival 
angle every predetermined angle step A9 using the evaluation function given by 
Eq.15. Here, the arrival-angle evaluation function given by Eq.14 is described in 
J. Capon, "High-Resolution Frequency- Wavenumber Spectrum Analysis." Proc. 
IEEE, 57 (8), pp. 1408 - 1418, 1969. 
15 Since elements outside the lower triangular part in lower triangular matrix 

L are null, the sum-of-product calculation amount for the arrival-angle evaluation 
function given by Eq.15 is ratio (M+3) / [2 (M+l)] lower than that for the 
conventional arrival-angle evaluation function given by Eq.14. A relation 
between the calculation amount for the cholesky factorization and that of Fc2(80 
20 for one arrival angle 9i given by Eq.15 is similar to that given by embodiment 1. 
Therefore, an increment of the calculation amount caused by the cholesky 
factorization can be considered sufficiently smaller than a decrement of the 
calculation amount for the arrival-angle evaluation in the entire angle range in 
Eq.15. 

25 Arrival-angle determination unit 23 detects a peak direction based on an 



13 

arrival-angle evaluation result every AG in a variable range of 6, and uses the 
detected direction as an arrival-angle estimation value of the incident waves. 

In the present embodiment, using arrival-angle evaluation function Fc2(8) 
(Eq.15), the calculation amount in the ai'iival- angle evaluation can be reduced 
5 compared with the arrival-angle evaluation function (Eq.14) by the Capon 
method. Function Fc2(0) has been derived by equivalent conversion of the 
arrival-angle evaluation function given by Eq.14, using lower triangular matrix L 
determined in triangular matrix calculation unit 21 by applying, the cholesky 
factorization to inverse matrix R- 1 of the correlation matrix. 

10 The direction estimation based on the Capon method has been discussed 

above. However, the method of the present embodiment can be as-is applied to a 
received signal sampled eveiy frequency when the delay time of the received 
wave is estimated at high resolution. That is because an evaluation function 
similar to that in Eq.14 is used. 

15 Additionally, coi'relation matrix calculation unit 7 can apply a spatial 

smoothing technique to the correlation matrix in order to suppress correlation 
wave. The method of the present embodiment can be similarly applied to this 
case, if the spatially smoothed correlation matrix instead of correlation matrix R 
is fed into the direction estimating processor. 

20 When an array antenna having the constant-interval linear array shape is 

employed, a steeling vector can be converted into a real vector using unitary 
matrix Qm given by Eq.12, thanks to the conjugation center symmetry of the 
phase of the steeling vector. In Eq.12, a(9) is a steering vector when the phase 
center matches with the array center. A method similar to the method in the 

25 present embodiment can be applied to this case, if direction estimating processor 
12a uses the real part of Qm h RQm instead of correlation matrix R, and b(8) 
instead of steering vector a(8). 
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When an array antenna having the linear array shape is employed, 
estimation accuracy in the bore-sight direction decreases. Therefore, arrival- 
angle evaluation unit 22 sets the angle interval in the end fire direction of the 
array antenna to be larger than that in the bore-sight direction, and calculates an 
5 evaluation value of the arrival-angle evaluation function. The calculation 
amount can be thus reduced. 

Additionally, triangular matrix calculation unit 21 may derive lower 
triangular matrix L and diagonal matrix D using the modified cholesky 
factorization. The modified cholesky factorization does not require the square 
10 root calculation, so that the calculation time can be reduced. 



3. THIRD EXEMPLARY EMBODIMENT 

Fig. 5 is a block diagram illustrating the other configuration of direction 
estimating px^ocessor 12 of the radio-wave arrival-direction estimating apparatus 
15 of the present invention. Components other than the direction estimating 
processor in the radio-wave arrival- dix^ection estimating apparatus are similar to 
those in embodiment 1 described in Fig. 1, so that diagrams and descriptions of 
these components are omitted. Components different from embodiment 1 will be 
mainly described hereinafter. 
20 Operations until correlation matrix R are fed into direction estimating 

processor 12b are similar to those in embodiment 1. 

Since correlation matrix R is a positive definite matrix, triangular matrix 
calculation unit 24 derives lower triangular matrix L given by 

R=LLH (Equation 16), 

25 using the cholesky factorization. 

Inverse matrix calculation unit 25 calculates inverse matrix Lr 1 of lower 
triangular matrix L. 
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Using lower triangular matrix L, arrival-angle evaluation unit 26 converts 
the arrival-angle evaluation function (Eq.14) derived by the Capon method to an 
evaluation function Fc3(0) expressed as 

Fc i(Q) - u — ?rr (Equation 17), 

3V J \\L- l a(6)\\ 2 1 

5 where, | |x| I is the norm of vector x, and a(9) is a normalized steering vector of 
the array antenna. Arrival-angle evaluation unit 26 then evaluates an arrival 
angle every predetermined angle step A9 using the evaluation function of Fca(9). 

Since elements outside the lower triangular part in lower triangular 
matrix L are null, the sum-of-product calculation amount for the arrival-angle 

10 evaluation function given by Eq.17 is ratio (M+3) / [2 (M+l)] lower than that for 
the Capon method's conventional arrival-angle evaluation function (Eq.14). A 
relation between the calculation amount (Eq.16) for the cholesky factorization 
and that of Fc2(0i) per arrival angle 8i is similar to that shown in embodiment 1. 
Therefore, an increment of the calculation amount caused by the cholesky 

15 factorization can be considered sufficiently smaller than a decrement of the 
calculation amount for the arrival-angle evaluation in the entire angle range in 
Eq.17. 

Arrival-angle determination unit 27 detects a peak direction based on an 
arrival-angle evaluation result every AG in a variable range of 9, and uses the 
20 detected direction as an arrival-angle estimation value of the incident waves. 

In the present embodiment, using arrival-angle evaluation function Fc3(9) 
(Eq.17), the calculation amount in the arrival-angle evaluation can be 
significantly reduced compared with the conventional arrival-angle evaluation 
function (Eq.14) by the Capon method. Function Fc3(9) has been derived by 
25 equivalent conversion of the arrival-angle evaluation function given by Eq.14, 
using lower triangular matrix L determined in triangular matrix calculation unit 
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24 by applying the cholesky factorization to correlation matrix R. 

The direction estimation based on the Capon method has been discussed 
above. However, the method of the present embodiment can be as-is applied to a 
received signal sampled every frequency when the delay time of the received 
5 wave is estimated at high resolution. That is because an evaluation function 
similar to that in Eq.14 is used. 

Additionally, correlation matrix calculation unit 7 can apply a spatial 
smoothing technique to the correlation matrix in order to suppress correlation 
wave. The method of the present embodiment can be similarly applied to this 

10 case, if the spatially smoothed correlation matrix instead of correlation matrix R 
is fed into the direction estimating processor. 

When an array antenna having the constant-interval linear array shape is 
employed, a steering vector can be converted into a real vector using unitary 
matrix Qm given by Eq.12, thanks to the conjugation center symmetry of the 

15 phase of the steering vector. In Eq.12, a(0) is a steering vector when the phase 
center matches with the array center. A method similar to the method in the 
present embodiment can be applied to this case, if direction estimating processor 
12b uses the real part of Qm h RQm instead of correlation matrix R, and b(6) 
instead of steering vector a(0). 

20 When an array antenna having the linear array shape is employed, 

estimation accuracy in the bore-sight direction decreases. Therefore, arrival- 
angle evaluation unit 26 sets the angle interval in the end fire direction of the 
array antenna to be larger than that in the bore-sight direction, and calculates an 
evaluation value of the arrival-angle evaluation function. The calculation 

25 amount can be thus reduced. 

Additionally, triangular matrix calculation unit 24, using the modified 
cholesky factorization, may derive lower triangular matrix L and diagonal matrix 
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D given by 

R = LD L H (Equation 18). 

The modified cholesky factorization does not require the square root calculation, 
so that the calculation time can be reduced. 

An arrival-angle evaluation function in this case is expressed as 



Pc 3 (fl ) = 



V' 



(Equation 19), 



where bk and dk are vector elements given by 

/ b x \ 

b 2 



L ~ 1 a (6 ) 



(Equation 20) 



and 



D = 



(d x 0 
0 d , 



0 



0 \ 
0 



(Equation 21), 



x-espectively. The lower triangular matrix has been used in the present 
embodiment; however, an upper triangular matrix also obviously produces a 
similar advantage. 



15 4. FOURTH EXEMPLARY EMBODIMENT 

Fig. 6 is a block diagram illustrating the other configuration of direction 
estimating processor 12 of the radio-wave arrival- direction estimating apparatus 
of the present invention. Components other than the direction estimating 
processor in the radio-wave arrival-direction estimating apparatus are similar to 

20 those in embodiment 1 described in Fig. 1, so that diagrams and descriptions of 
these components are omitted. Components different from embodiment 1 will be 
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mainly described hereinafter. 

Operations until correlation matrix R are fed into direction estimating 
processor 12c is similar to those in embodiment 1. 

Since correlation matrix R is a positive definite matrix, triangular matrix 
calculation unit 28, using the cholesky factorization, derives lower triangular 
matrix L given by 

R = LL« (Equation 22). 

Using lower triangular matrix L, arrival-angle evaluation unit 29 converts 
conventional arrival-angle evaluation function Ff(6) obtained by Fourier method 
as given by 

F F {6) = a H (8)Ra(6) (Equation 23), 

to an evaluation function expressed as 

F F2&) HI L H a{6) || 2 (Equation 24), 

where, I I x | | is the norm of vector x, and a(8) is a normalized steering vector of 
the array antenna. Arrival-angle evaluation unit 29 then evaluates an arrival 
angle every predetermined angle step A6 using the evaluation function given by 
Eq.24. Here, the arrival- angle evaluation function given by Eq.23 is described in 
M. S. Bartlett. "Smoothing Periodograms from Time Series with Continuous 
Spectra." Nature, 161, pp. 686 - 687, (1948). 

Since elements outside the lower triangular part in lower triangular matrix 
L are null, the sum-of-product calculation amount for the arrival-angle evaluation 
function given by Eq.24 is ratio (M+3) / [2 (M+l)] lower than that for the 
conventional arrival-angle evaluation function given by Eq.23. An increment of 
the calculation amount caused by the cholesky factorization is sufficiently smaller 
than a decrement of the calculation amount for the arrival-angle evaluation in the 
entire angle range in Eq.24. That is because the relation between the 
calculation amounts in the present embodiment is also similar to that shown in 
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embodiment 1. 

Arrival-angle determination unit 30 detects a peak direction based on an 
arrival-angle evaluation result eveiy A6 in a variable range of 6, and uses the 
detected direction as an arrival-angle estimation value of the incident waves. 

In the present embodiment, using arrival-angle evaluation function Ff2(6) 
(Eq.24), the calculation amount in the arrival-angle evaluation can be 
significantly reduced compared with the conventional arrival-angle evaluation 
function (Eq.23) by Fourier method. Function Ff2(8) has been derived by 
equivalent conversion of the arrival-angle evaluation function given by Eq.23, 
using lower triangular matrix L determined in the triangular matrix calculation 
unit by applying the cholesky factorization to the correlation matrix. 

Additionally, correlation matrix calculation unit 7 can apply a spatial 
smoothing technique to the correlation matrix in order to suppress correlation 
wave. The method of the present embodiment can be similarly applied to this 
case, if the spatially smoothed correlation matrix instead of correlation matrix R 
is fed into direction estimating processor 12c. 

When an array antenna having the constant-interval linear array shape is 
employed, a steering vector can be converted into a real vector using unitary 
matrix Qm given by Eq.12, thanks to the conjugation center symmetry of the 
phase of the steering vector. In Eq.12, a(0) is a steering vector when the phase 
center matches with the array center. A method similar to the method in the 
present embodiment can be applied to this case, if direction estimating processor 
12c uses the real part of Qm h RQm instead of correlation matrix R, and b(0) 
instead of steering vector a(6). 

When an array antenna having the linear array shape is employed, 
estimation accuracy in the bore-sight direction decreases. Therefore, arrival- 
angle evaluation unit 30 sets the angle interval in the end fire direction of the 
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array antenna to be larger than that in the bore-sight direction, and calculates an 
evaluation value of the arrival-angle evaluation function. The calculation 
amount can be thus reduced. 

Additionally, triangular matrix calculation unit 28, using the modified 
cholesky factorization, may derive lower triangular matrix L and diagonal matrix 
D given by Eq.18. The modified cholesky factorization does not require the 
square root calculation, so that the calculation time can be reduced. An 
arrival-angle evaluation function in this case is expressed as 



M u2 

where bk and dk are vector elements written as 
= L"a{d) 



(Equation 25), 



(Equation 26) 



and 



Id, 



D 



0 

d 2 



0 



0 



0 \ 
0 

' M J 



(Equation 27), 



respectively. The lower triangular matrix has been used in the present 
embodiment; however, an upper triangular matrix also obviously produces a 
similar advantage. 



5. FIFTH EXEMPLARY EMBODIMENT 

Fig. 7 is a block diagram illustrating the other configuration of direction 
20 estimating processor 12 of the radio-wave arrival-direction estimating apparatus 
of the present invention. Components other than the direction estimating 
processor in the radio-wave arrival-direction estimating apparatus are similar to 
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those in embodiment 1 described in Fig. 1, so that diagrams and descriptions of 
these components are omitted. Components different from embodiment 1 will be 
mainly described hereinafter. 

Operations until complex digital signals 6 are obtained are similar to those 
in embodiment 1. 

Correlation vector calculation unit 31 selects one of complex digital signals 
6, and performs a correlation calculation between an antenna element — a 
reference antenna — corresponding to the selected signal and another antenna 
element, thereby deriving a correlation vector. An example will be described 
hereinafter employing antenna element 1-1 as the reference antenna. Antenna 
elements 1-1 to 1-M receive sampling data xi (k), X2 (k), ... , xm (k) at sampling 
time to + kAT, respectively. Here to is an arbitrary time, AT is a sampling 
interval, and k is a natural number. Correlation vector calculation unit 31, 
using snapshot data of the sampling data for k = 1 to N, derives M-dimensional 
correlation vector Rv expressed as 



where asterisk * shows complex conjugate. 

Arrival-angle evaluation unit 32 evaluates an arrival angle every 
predetermined angle step A0 using the evaluation function given by 



where, I I x | I is the norm of vector x, and a(0) is a normalized steering vector of 
the array antenna. 

Arrival-angle determination unit 33 detects a peak direction based on an 
arrival-angle evaluation result every A9 in a variable range of 9, and uses the 
detected direction as an arrival-angle estimation value of the incident waves. 




(Equation 28), 



F v (d) = \\R v H a(6)\\ 



(Equation 29), 
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The arrival angle is evaluated using the correlation vector instead of the 
correlation matrix in the present embodiment, so that the calculation amount in 
the arrival-angle evaluation can be significantly reduced compared with the 
conventional Fourier method (Eq.23). 
5 When an array antenna having the linear array shape is employed, 

estimation accuracy in the bore-sight direction decreases. Therefore, arrival- 
angle evaluation unit 32 sets the angle interval in the end fire direction of the 
array antenna to be larger than that in the bore-sight direction, and calculates an 
evaluation value of the arrival-angle evaluation function. The calculation 
10 amount can be thus reduced. 

6. SIXTH EXEMPLARY EMBODIMENT 

Fig. 8 is a block diagram illustrating the other configuration of direction 
estimating processor 12 of the radio-wave arrival-direction estimating apparatus 

15 of the present invention. Components other than the direction estimating 
processor in the radio-wave arrival-direction estimating apparatus are similar to 
those in embodiment 1 described in Fig. 1, so that diagrams and descriptions of 
these components are omitted. Components different from embodiment 1 will be 
mainly described hereinafter. The present embodiment employs an array 

20 antenna having the constant-interval linear array shape. 

Operations until complex digital signals 6 are obtained are similar to those 
in embodiment 1. Correlation vector calculation unit 34 selects one of complex 
digital signals 6, and performs a correlation calculation between an antenna 
element ~ a reference antenna — corresponding to the selected signal and another 

25 antenna element, thereby deriving a correlation vector. An example will be 
described hereinafter employing antenna element 1-1 as the reference antenna. 
Antenna elements 1-1 to 1-M receive sampling data xi (k), X2 (k), ... , xm (k) at 
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sampling time to + kAT, respectively. Here to is an arbitraiy time, AT is a 
sampling interval, and k is a natural number. Correlation vector calculation 
unit 34, using snapshot data of the sampling data for k = 1 to N, derives M- 
dimensional correlation vector Rv expressed as Eq.28. In Eq.28, asterisk * 
shows complex conjugate. 

Since the array antenna has the constant-interval linear array shape, a 
steering vector can be converted into a real vector using unitary matrix Qm given 
by Eq.12, thanks to the conjugation center symmetiy of the phase of the steering 
vector. In Eq.12, a(0) is a steering vector when the phase center matches with 
the array center. Unitary transforming unit 35 unitary-transforms correlation 
vector Rv as in 

R r ! =real (q* M R vQ M ) (Equation 30), 

R r 2 =real(q*i, m R v Q Nt ) (Equation 31), 

where qy is (i, j) element of matrix Qm, real (x) is a vector comprising real parts of 
respective elements of vector x, and m is M / 2 + 1 for even M (number of 
elements), or m is (M + 1) / 2 + 1 for odd M. 

Arrival-angle evaluation unit 36 evaluates an arrival angle every 
predetermined angle step A8 using the evaluation function given by 

F*2&) - [bx(0)R rl + b m (6)R r2 ]b(d) (Equation 32), 

where, real steering vector b(0) is converted from steering vector a(9) using Eq.12, 
and bk(0) is the k-th element of real steering vector b(0). 

Arrival-angle determination unit 37 detects a peak direction based on an 
arrival-angle evaluation result every A0 in a variable range of 0, and uses the 
detected direction as an arrival-angle estimation value of the incident waves. 

The arrival angle is evaluated using the correlation vector instead of the 
correlation matrix and using the real steering vector in the present embodiment, 
so that the calculation amount in the arrival-angle evaluation can be significantly 
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reduced compared with the conventional Fourier method (Eq.23). 

When an array antenna having the linear array shape is employed, 
estimation accuracy in the bore-sight direction decreases. Therefore, arrival- 
angle evaluation unit 37 sets the angle interval in the end fire direction of the 
5 array antenna to be larger than that in the bore-sight direction, and calculates an 
evaluation value of the arrival-angle evaluation function. The calculation 
amount can be thus reduced. 

7. SEVENTH EXEMPLARY EMBODIMENT 
10 Fig. 9 is a block diagram illustrating the other configuration of arrival- 

angle evaluation unit 10 of the radio-wave arrival-direction estimating apparatus 
of the present invention. Components other than the arrival-angle evaluation 
unit in the radio-wave arrival-direction estimating apparatus are similar to those 
in embodiment 1 described in Fig. 1, so that diagrams and descriptions of these 
15 components are omitted. Components different from embodiment 1 will be 
mainly described hereinafter. 

The present embodiment assumes that antenna elements 1-1 to 1-M of 
array antenna 1 are arranged linearly at a constant interval. Operations until 
triangular matrix L is fed into arrival-angle evaluation unit 10a are similar to 
20 those in embodiment 1. 

Arrival-angle evaluation unit 10a comprises the following components: 
positive-region evaluation unit 40 for calculating an arrival-angle 
evaluation function in the positive angle region (0° ^ 0 s 90°) using the bore-sight 
direction of the array antenna as an angle reference (9 = 0), because the array 
25 antenna has the constant-interval linear array shape; and 

negative-region evaluation unit 41 for converting the evaluation result of 
the positive-region evaluation unit to an arrival-angle evaluation value in the 
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negative angle region (-90° <; 6 <s 0°). 

Since the array antenna has the constant-interval Linear array shape, steering 

vector a(8) is a complex vector expressed as 
exp{- jljzd ■ 0 * sin 6 I x} 



a(8) = 



exp{- j2itd • 1 - sin 6 I x} 



(Equation 33). 



exp{- jlnd - (M - 1) • sin 6 I x) 

5 The real part of the complex vector is an even function of 0, and the imaginary 

part is an odd function of 6. Using this feature, Eq.7 can be transformed to 

F 2 (±0,)4^V^)li]" 2 

- [| (Rc(L H )+ ; lm{L H ))- (Re( fl (0,))± ; lm(a(0,)))||} 2 

= [|| Re{L" )• Re( 0 (6 f ))T Imfc" > lm(a(0,))|| 2 - || Imfc" )• Re(n(e,))± Re(l" )- Im^C^ ))l| 2 ]~ J 
= [||clTc2|| 2 -||c3^c4|rj- 1 

(Equation 34), 

for 6i satisfying 0° <; 8i =s 90°. Here, Re(x)'is a vector comprising real parts of 
10 respective elements of the complex vector (or matrix) x, Im(x) is a vector 
comprising imaginary parts of them, d is an interval of antenna elements, K is a 
wavelength of carrier frequency, and vectors cl, c2, c3, c4 are given by 

cl = Re(z," )• Re(a(0,)) (Equation 35), 

c2 = Im^J-ImCflCfl,)) (Equation 36), 

15 c3 = lm(L H )• Re(a(0,)) (Equation 37), 

and 

c4 - Re(L" )• lm(a(0, )) (Equation 38), 

respectively. 

Positive-region evaluation unit 40 calculates an arrival-angle evaluation 
20 function in the positive angle region (0° 6 =s 90°) at predetermined angle step A0 
using the relation discussed above. Positive-region evaluation unit 40 firstly 
calculates vectors cl, c2, c3, c4, and then derives a positive-region arrival-angle 
evaluation value in accordance with 
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F3(9,) = [|| cl - c2 || 2 - || c3 + c4 || 2 J* 1 (Equation 39). 

Evaluation unit 40 outputs resultant values of vectors cl, c2, c3, c4 to negative- 
region evaluation converter 41. 

Negative-region evaluation converter 41 calculates an arrival-angle 
5 evaluation value in the -6i direction using vectors cl, c2, c3, c4 supplied from 
evaluation unit 40 in accordance with 

^ 4 (-S/) - [II cl + c2 || 2 - || c3 - c4 If 8 ]" 1 (Equation 40). 

Arrival-angle determination unit 11 detects a peak direction based on an 
arrival-angle evaluation result of arrival-angle evaluation unit 10a every AG in a 
10 variable range of 6, and uses the detected direction as an arrival-angle estimation 
value of the incident waves. 

In the present embodiment, the positive -region arrival-angle evaluation 
value can be converted to the negative-region arrival-angle evaluation value 
using vectors cl, c2, c3, c4, when the array antenna having the constant : interval 
15 linear array shape is employed. Here the vectors cl, c2, c3, c4 are calculated 
when the positive -region arrival-angle evaluation value is derived. The 
calculation amount in the arrival-angle evaluation can be further reduced 
substantially in half. 

The present embodiment has been described using the arrival-angle 
20 evaluation function given by Eq.7 in embodiment 1. When an array antenna 
having the constant-interval linear array shape is employed, a similar 
transformation can be also applied to the arrival-angle evaluation function using 
steering vector a(0) in the other embodiment. The conversion to the negative- 
region arrival-angle evaluation value is allowed using vectors cl, c2, c3, c4 
25 calculated when the positive-region arrival-angle evaluation value are derived, as 
discussed above. The calculation amount in the arrival-angle evaluation can be 
reduced substantially in half compared with the conventional method. 
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8. EIGHTH EXEMPLARY EMBODIMENT 

Fig. 10 is a block diagram illustrating the other configuration of direction 
estimating processor 12 of the radio-wave arrival- direction estimating apparatus 
5 of the present invention. Components other than the direction estimating 
processor in the radio-wave arrival-direction estimating apparatus are similar to 
those in embodiment 1 described in Fig. 1, so that diagrams and descriptions of 
these components are omitted. Components different from embodiment 1 will be 
mainly described hereinafter. 

10 The present embodiment assumes that antenna elements 1-1 to 1-M of 

array antenna 1 are arranged linearly at a constant interval. Operations until 
correlation matrix calculation unit 7 outputs correlation matrix R are similar to 
those in embodiment 1. 

Unitary transforming unit 42 transforms correlation matrix R as given by 

15 R r =real(Q M H R Qm) (Equation 41), 

using unitary matrix Qm for converting a steering vector employing the phase 
center as the array center to a real vector. Hex*e real(A) is a matrix comprising 
real parts of respective elements of matrix A. 

Noise spatial eigenmatrix calculation unit 8a applies the eigenvalue 

20 factorization to unitary-transformed correlation matrix Rr by unitary 
transforming unit 42 to derive eigenvalues Xi - Xu in the descending order and 
eigenvactors ei - eM corresponding to them. When a number of the incident 
waves is S, calculation unit 8a outputs noise spatial eigenmatrix En = [es+i, , 
eM] comprising (M - S) eigenvactors, namely column vectors, belonging to a noise 

25 partial space having the relation given by Eq.3. 

Triangular matrix calculation unit 9a derives product U of a noise spatial 
eigenmatrix and a complex-conjugated and transposed matrix of it as in Eq.5. 
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Triangular matrix calculation unit 9a derives lower triangular matrix L given by 
Eq.6 using the cholesky factorization, because matrix U is an M x M positive 
definite matrix. 

Arrival-angle evaluation unit 10a comprises the following components: 
positive-region evaluation unit 40a for calculating the arrival-angle 

evaluation function in the positive angle region (0° ^ 0 ^ 90°) using the bore-sight 

direction of the array antenna as an angle reference (9 = 0), because the array 

antenna has the constant-interval linear array shape; and 

negative-region evaluation unit 41a for converting the evaluation result 

of the positive -region evaluation unit to an arrival-angle evaluation value in the 

negative angle region (-90° 0°). 

When steering vector a(8) employing the phase center as the array center is 
converted using unitary matrix Qm, real vector b(0) is derived as in 



b(6) = V2 



cos^ M 2 1 f.i j, • ■ • , cos(/,< ), - sin ^ M 2 1 A' j, * • ' , - sin(,a ) 



for M = 2m 



or 



b(6) = V2 



for M = 



2m 4- 1 , (Equation 42), 

where \i is written as 

\jl = - ^-d sin# (Equation 43), 

A 

where, d is an interval between the antenna elements, and X is a wavelength of 
carrier frequency. 

When number M of antenna elements equals 2m, b(6) is an even function 
for elements 1 to m or an odd function for elements m+1 to 2m, as shown in Eq.42. 
Using this feature, Eq.7 can be transformed to 
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15 



F2( ± e,)4\L"b(±e l )\\\ 1 

4^(0±^(*,)ll]" 2 
= (j|cl±c2||]- 2 



(Equation 44), 



for 9i satisfying 0° s. 6i s 90°. Here, cl, c2, b eve n(e), and b 0 dd(9) are given by 
cl-^^Cc,) (Equation 45), 

cl-H'b^ie,) (Equation 46), 



^ n (0) = V2 



COS 



( M -1 



/.i^-scos^OA-sO 



(Equation 47), 



and 



/ Af - 1 \ / x 

0,---,0,-sin ,u ,* * • 7 -sin( j aJ 

V 2 / 



(Equation 48), 



respectively. 

When number M of antenna elements equals 2m + 1, b(6) is an even 
10 function for elements 1 to m, and b(0) is an odd function for elements m+2 to M. 
Using this feature, Eq.7 can be transformed to Eq.44, for 6i satisfying 0° <; 8i ^ 90°. 
In this case, b e ven(9) and b 0 dd(8) are given by 

T 



^„(0) = V2 



and 



0, ■ • • ,0,- sin^ ^^ j, • • • sin(jii ) 



(Equation 49), 



(Equation 50), 



respectively. 

Positive-region evaluation unit 40a calculates an arrival-angle evaluation 
function in the positive angle region (0° <; G <; 90°) at predetermined angle step AO 
using the relation discussed above. Positive -region evaluation unit 40a firstly 
20 calculates vectors cl, c2 (Eq.45 and Eq.46) and then derives a positive-region 
arrival-angle evaluation value in accordance with 
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F3(0,) = [|| cl + c2 II 2 ]" 1 (Equation 51). 

Evaluation unit 40a outputs resultant values of vectors cl, c2 to negative-region 
evaluation converter 41a. 

Negative-region evaluation converter 41a calculates an arrival-angle 
5 evaluation value in the -8i direction using vectors cl, c2 supplied from evaluation 
unit 40a in accordance with 

F4(-0 / ) = [|| cl - c2 I) 2 ]" 1 (Equation 52). 

Arrival-angle determination unit 11a detects a peak direction based on an 
arrival-angle evaluation result of arrival-angle evaluation unit 10a every A6 in 
10 the variable range of 0, and uses the detected direction as an arrival-angle 
estimation value of the incident waves. 

In the present embodiment, the array antenna having the constant- 
interval linear array shape has been employed, and the real steering vector has 
been derived by unitary matrix transformation. The positive-region arrival- 
15 angle evaluation value can be converted to the negative-region arrival-angle 
evaluation value, using vectors cl, c2 calculated when the positive -region arrival- 
angle evaluation value are derived. The calculation amount in the arrival-angle 
evaluation can be reduced substantially in half compared with embodiment 1. 

20 9. NINETH EXEMPLARY EMBODIMENT 

Fig. 11 is a block diagram illustrating the other configuration of direction 
estimating processor 12 of the radio-wave arrival- direction estimating apparatus 
of the present invention. Components other than the direction estimating 
processor in the radio-wave arrival-direction estimating apparatus are similar to 

25 those in embodiment 1 described in Fig. 1, so that diagrams and descriptions of 
these components are omitted. Components different from embodiment 1 will be 
mainly described hereinafter. 
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The present embodiment assumes that antenna elements 1-1 to 1-M of 
array antenna 1 are arranged linearly at a constant interval. Operations until 
correlation matrix calculation unit 7 outputs correlation matrix R are similar to 
those in embodiment 1. 
5 Unitary transforming unit 42 transforms correlation matrix R as given by 

Eq.41, using unitary matrix Qm for converting a steering vector employing the 
phase center as the array center to a real vector. Inverse matrix calculation unit 
20a calculates an inverse matrix of unitary-transformed correlation matrix Rr 
derived by the unitary transforming unit. Triangular matrix calculation unit 
10 21a derives lower triangular matrix L given by Eq.13 using the cholesky 
factorization. That is becaxise inverse matrix R 1 is a positive definite matrix. 
Arrival-angle evaluation unit 22a comprises the following components: 
positive-region evaluation unit 53 for calculating an arrival-angle 
evaluation function in the positive angle region (0° <. 9 <; 90°) using the bore-sight 
15 direction of the array antenna as an angle reference (0 = 0), because the array 
antenna has the constant-interval linear array shape; and 

negative -region evaluation unit 54 for converting the evaluation result of 
the positive -region evaluation unit to an arrival- angle evaluation value in the 
negative angle region (-90° <s 0 <; 0°). 
20 Positive-region evaluation unit 53 calculates an arrival-angle evaluation function 
in the positive angle region (0° <; 0 ^ 90°) at predetermined angle step A0. 
Positive-region evaluation unit 53 firstly calculates vectors cl, c2 given by Eq.45 
and Eq.46 for 0i satisfying 0° <> 0i <, 90° and then derives a positive-region 
arrival-angle evaluation value in accordance with 

25 

F ci( 0 t) - [II cl + c2 IP]" 1 (Equation 53). 

Evaluation unit 53 outputs resultant values of vectors cl, c2 to negative-region 
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evaluation converter 54. 

Negative-region evaluation converter 54 calculates an arrival-angle 
evaluation value in the -0i direction using vectors cl, c2 supplied from evaluation 
unit 53 in accordance with 

F c2(- 6 i) = [|| cl - c2 I) 2 ]" 1 (Equation 54). 

Arrival-angle determination unit 23 detects a peak direction based on an 
arrival-angle evaluation result of arrival-angle evaluation unit 22a eveiy A9 in 
the variable range of 0, and uses the detected direction as an arrival-angle 
estimation value of the incident waves. 

In the present embodiment, the array antenna having the constant- 
interval linear array shape has been employed, and the real steering vector has 
been derived by unitary matrix transformation. The positive-region arrival- 
angle evaluation value can be converted to the negative -region arrival-angle 
evaluation value using vectors cl, c2 calculated when the positive-region arrival- 
angle evaluation value are derived. The calculation amount in the arrival-angle 
evaluation can be reduced substantially in half compared with embodiment 1. 

The present embodiment can be similarly applied to the arrival-angle 
evaluation function described in embodiment 3. 

10. TENTH EXEMPLARY EMBODIMENT 

Fig. 12 is a block diagram illustrating the other configuration of direction 
estimating processor 12 of the radio-wave arrival- direction estimating apparatus 
of the present invention. Components other than the direction estimating 
processor in the radio-wave arrival- direction estimating apparatus are similar to 
those in embodiment 1 described in Fig. 1, so that diagrams and descriptions of 
these components are omitted. Components different from embodiment 1 will be 
mainly described hereinafter. 
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The present embodiment assumes that antenna elements 1-1 to 1-M of 
array antenna 1 are arranged linearly at a constant interval. Operations until 
correlation matrix calculation unit 7 outputs correlation matrix R are similar to 
those in embodiment 1. 
5 Unitary transforming unit 42 transforms correlation matrix R as given by 

Eq.41, using unitary matrix Qm for converting a steering vector employing the 
phase center as the array center to a real vector. Triangular matrix calculation 
unit 24a derives lower triangular matrix L using the cholesky factorization, 
because unitary-transformed correlation matrix Rr is a positive definite matrix. 

10 Arrival- angle evaluation unit 25a comprises the following components: 

positive -region evaluation unit 55 for calculating an arrival-angle 
evaluation function in the positive angle region (0° <; 0 ' <, 90°) using the bore-sight 
direction of the array antenna as an angle reference (6 = 0), because the array 
antenna has the linear array shape; and 

15 negative-region evaluation unit 56 for converting the evaluation result of 

the positive-region evaluation unit to an arrival-angle evaluation value in the 
negative angle region (-90° «s 6 <; 0°). 

Positive-region evaluation unit 55 calculates an arrival-angle evaluation function 

in the positive angle region (0° <; 0 =s 90°) at predetermined angle step A0. 
20 Positive -region evaluation unit 55 firstly calculates vectors cl, c2 given by Eq.45 

and Eq.46 for 0i satisfying 0° ^ 6i ^ 90° and then derives a positive -region 

arrival-angle evaluation value in accordance with 

- II cl + c2 IP (Equation 55). 

Evaluation unit 55 outputs respective resultant values of vectors cl, c2 to 
25 negative-region evaluation convei'ter 56. 

Negative-region evaluation converter 56 calculates an arrival-angle 

evaluation value in the -8i direction using vectors cl, c2 supplied from evaluation 



34 



unit 55 in accordance with 

F F2 (-6 t ) = || cl - c2 || 2 (Equation 56). 

Arrival-angle determination unit 26a detects a peak direction based on an 
arrival-angle evaluation result of arrival-angle evaluation unit 25a every A9 in 
5 the variable range of 6, and uses the detected direction as an arrival-angle 
estimation value of the incident waves. 

In the present embodiment, the array antenna having the constant- 
interval linear array shape has been employed, and the real steering vector has 
been derived by unitary matrix transformation. The positive-region arrival- 
10 angle evaluation value can be converted to the negative-region arrival-angle 
evaluation value using vectors cl, c2 calculated when the positive-region arrival- 
angle evaluation value are derived. The calculation amount in the arrival-angle 
evaluation can be reduced substantially in half compared with embodiment 1. 

The present embodiment can be applied to embodiment 6 by a similar 
15 equation transformation. 

11. ELEVENTH EXEMPLARY EMBODIMENT 

Fig. 13 is a block diagram illustrating the other config-uration of direction 

estimating processor 12 of the radio-wave arrival-direction estimating apparatus 
20 of the present invention. Components other than the direction estimating 

processor in the radio-wave arrival-direction estimating apparatus are similar to 

those in embodiment 1 described in Fig. 1, so that diagrams and descriptions of 

these components are omitted. Components different from embodiment 1 will be 

mainly described hereinafter. 
25 The configuration in the present invention includes the following 

components, in addition to direction estimating processor 12 shown in Fig. 1: 

high-accuracy arrival-angle evaluation unit 60 for calculating an 
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evaluation value of an arrival-angle evaluation function at an angle interval 
smaller than an angle interval calculated by an arrival-angle evaluation unit; and 
high-accuracy arrival-angle determination unit 61 for highly accurately 
determining an arrival angle based on the evaluation value by high-accuracy 
arrival-angle evaluation unit 60. 

The following operations are similar to those in embodiment 1: 

correlation matrix R is fed into direction estimating processor 12k; 
triangular matrix calculation unit 9 calculates triangular matrix L; 

and 

arrival-angle determination unit 11 detects a peak direction based on an 
arrival-angle evaluation result of arrival- angle evaluation unit 10 every A8 in a 
variable range of 6, and uses the detected direction as an arrival-angle estimation 
value of the incident waves. 

Fig. 14 is a graph illustrating operations of high-accuracy arrival-angle 
evaluation unit 60 and high-accuracy arrival-angle determination unit 61. 
Arrival-angle evaluation unit 60 reevaluates k arrival angles 0k (k is a natural 
number) supplied from arrival-angle determination unit 11 in angle range Wk of 
<j) satisfying (0k - A0) < <j) < (0 + A0), at angle step A(j> smaller than the angle step AO 
in arrival-angle evaluation unit 10, and using the arrival-angle evaluation 
function. Arrival-angle evaluation unit 60 then outputs the resultant arrival 
angle, which is highly accurate, to arrival-angle determination unit 61. 

Arrival-angle determination unit 61, based on evaluated values in k angle 
ranges Wk, detects peak direction <|)peak,k in each angle range Wk, and outputs 
the peak direction as an evaluated high-accuracy arrival angle. 

In the present embodiment, an arrival-angle has been more accurately 
reevaluated at a step smaller than angle step A0 in arrival-angle evaluation unit 
10 restrictively in a range around the arrival angle estimated in arrival-angle 
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determination unit 11. Therefore, the direction estimating processor in the 
present embodiment not only produces the advantage shown in embodiment 1, 
but also can highly accurately estimate the arrival angle without unnecessarily 
increasing the calculation amount for the arrival-angle evaluation function. 
5 The operations in the present embodiment have been described in the 

structure comprising high accuracy evaluation unit 60 and high-accuracy arrival- 
angle determination unit 61 in addition to direction estimating processor 12 
shown in embodiment 1. The present embodiment can produce a similar 
advantage even in a structure comprising these units in addition to direction 
10 estimating processor shown in each of embodiments 2 to 10. 

12. TWELFTH EXEMPLARY EMBODIMENT 

Fig. 15 is a block diagram illustrating a structure of a directivity variable 
receiver of the present invention. The directivity variable receiver in Fig. 15 

15 selects a plurality of sector antennas with different main beam directions, and 
changes the directivity, based on arrival-direction estimation result 64 derived by 
arrival-direction estimating apparatus 63 as described in embodiments 1 to 11. 
Operations in arrival-direction estimating apparatus 63 are similar to those 
described above, so that descriptions of the operations are omitted. Additional 

20 components will be described hereinafter. 

The directivity variable receiver comprises arrival-direction estimating 
apparatus 63, m 0>2) sector antennas 65-1 to 65-m with different main beam 
directions, sector switch 66, sector control unit 67, and receiving unit 68. 

Operations for estimating a radio-wave arrival direction using received 

25 signals 2-1 to 2-M obtained by an array antenna are similar to those described in 
embodiments 1 to 11. Here the array antenna comprises a plurality of antenna 
elements 1-1 to 1-M. Final arrival-direction estimation result 64 is fed into 
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sector control unit 67. Sector control unit 67, based on estimation result 64, 
selects the ms-th sector antenna having a main beam direction closest to the 
estimated direction from the plurality of sector antennas 65-1 to 65-m. Sector 
control unit 67 further controls sector switch 66 based on sector control signal 69 
5 to connect sector switch 66 to receiving unit 68. Receiving unit 68 demodulates 
signal 70 received by the ms-th sector antenna. 

These operations allow the selection of the optimal sector antenna having 
the main beam direction closest to the arrival direction from a plurality of. sector 
antennas 65-1 to 65-m, and provide received signal 70 with a high signal-to-noise 
10 ratio. Many waves having directions different from the main beam direction of 
the selected antenna can be suppressed, and interference between codes can be 
reduced. 

The present embodiment shows a structure of the receiver; however, the 
receiver can be also used as a transmitter by replacing receiving unit 68 with a 
15 transmitting unit. In this case, transmitted power is reduced and radio waves 
are not radiated in unnecessary directions, so that interference with other station 
can be reduced. 

Receiving unit 68 and transmitting unit 71 may be inter-coupled through a 
switch 72 as shown in Fig. 16 to provide a transceiver capable of 

20 switching between the reception and the transmission. 

13. THIRTEENTH EXEMPLARY EMBODIMENT 

Fig. 17 is a block diagram illustrating another structure of a du*ectivity 
variable receiver of the present invention. The directivity variable receiver in 
25 Fig. 17 changes directivity and receives signals, based on arrival-direction 
estimation result 64 derived by arrival-direction estimating apparatus 63 as 
described in embodiments 1 to 11. Operations in arrival-direction estimating 
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apparatus 63 are similar to those described in embodiments 1 to 11, so that 
descriptions of the operations are omitted. Additional components will be 
described hereinafter. 

L (>1) antenna elements 75-1 to 75-L may have a structure for dividing 
signals sent from antenna elements 1-1 to 1-M m Fig. 1 into two, but antenna 
elements having a different structure will be described hereinafter. Receiving 
units 77-1 to 77-L connected to respective antenna elements 75-1 to 75-L convert 
freqtiencies of signals 76-1 to 76-L received by antenna elements 75-1 to 75-L, and 
then demodulate the converted signals to signals 78-1 to 78-L comprising 
orthogonal I and Q signals. A/D converters 79-1 to 79-L convert respective 
demodulated signals 78-1 to 78-L, which are analog, to complex digital signals 
80-1 to 80-L. Sampling frequencies fs of A/D converters 79-1 to 79-L must satisfy 
Nyquist condition, fs > 2WB, in band WB (Hz) of transmitted modvriated wave so 
that the signals can be subsequently demodulated. Directivity control unit 81 
assigns complex weights to complex digital signals 80-1 to 80-L, based on arrival- 
direction estimation result 64 from arrival-direction estimating apparatus 63. 
Here the complex weights are vised for controlling the directivity. Receiving unit 
68 receives the weighted signals. 

As discussed above, directivity control unit 81 assigns complex weights to 
complex digital signals 80-1 to 80-L and combines the signals with each other to 
generate directivity in the direction of arrival-direction estimation result 64. A 
more optimal directivity pattern can be created in this case, communication 
qualit}' higher than that using a sector antenna is allowed. 

The present embodiment shows a structure for controlling the directivity in 
the receiver; however, a structure for controlling the. directivity in a transmitter 
may be also employed as shown in Fig. 18. In this case, transmitted power is 
reduced and radio waves are not radiated in unnecessary directions, so that 
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interference with other station can be reduced. 

In Fig. 18, directivity control unit 83 receives signal 82 transmitted from 
transmitting unit 71. Directivity control unit 83 divides transmitted signal 82 
into L signals. Directivity control unit 83 then assigns complex weights for 
5 controlling the directivity to respective divided signals 82, based on arrival- 
direction estimation result 64 from arrival-direction estimating apparatus 63, 
and outputs resultant complex digital signals 84-1 to 84-L. D/A converters 85-1 
to 85-L convert the digital signals to analog signals, and output the analog signals 
as base band signals 86-1 to 86-L. Transmission frequency converters 87-1 to 

10 87-L convert frequencies of base band signals 86-1 to 86-L to a transmission 
frequency band, and output resultant RF signals 88-1 to 88-L. 

Antenna elements 89-1 to 89-L transmit signals. 

A transceiver having functions shown in Fig. 17 and Fig. 18 can be 
employed. In this case, communication quality can be improved, and 

15 transmitted power is reduced and radio waves are not radiated in unnecessary 
directions to reduce interference with other station. 

The radio-wave arrival-direction estimating apparatus of the present 
invention can reduce a total calculation amount for arrival-angle evaluation using 
an arrival-angle evaluation function, without causing estimation accuracy 

20 degradation of the arrival direction. The estimating apparatus can also speed a 
calculation process or simplify an apparatus structure. Additionally, high 
quality communication is allowed, when the transmitting units or the receiving 
units in the transceiver have additional antennas having a directivity control 
function for generating directivity to the arriving direction of the arrival-direction 

25 estimating apparatus. 



